Detection of avoided crossings by fidelity 
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The fidelity, defined as overlap of eigenstates of two slightly different Hamiltonians, is proposed 
as an efficient detector of avoided crossings in the energy spectrum. This new application of fidelity 
is motivated for model systems, and its value for analyzing complex quantum spectra is underlined 
by applying it to a random matrix model and a tilted Bose-Hubbard system. 
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Introduction. - The enormous progress in cooling and 
manipulating ultracold atomic gases in recent years has 
opened new perspectives on interacting many-body mod- 
els from condensed matter physics Qj. It led to ques- 
tions and opportunities beyond conventional solid-state 
physics, e.g., the direct experimental study of quantum 
phase transitions [l[ , the role and engineering of genuine 
quantum correlations in such systems [H, and the phe- 
nomenon of quantum chaos in systems that consist of in- 
distinguishable particles which lack an obvious classical 
counterpart [3401 • In this context, it is possible to de- 
tect a quantum phase transition by the change of fidelity 
(modulus of the overlap between eigenstates of slightly 
different Hamiltonians) [tJ, since the ground state of a 
quantum system changes dramatically at a critical pa- 
rameter 0. 

Up to now, the temporal change of fidelity - as the over- 
lap of the same initial states evolved by different Hamil- 
tonians [9fl — has been measured experimentally in wave 
billiards [l(| , but also in systems of cold atoms subject to 
optical potentials [ill, E3 ■ Similar techniques may be ap- 
plied to measure the evolving overlap of two eigenstates 
where time is substituted by the change of some tunable 
control parameter. Often a quantum phase transition 
may be viewed, for finite-size realizations of a system, as 
an avoided crossing (AC) in parameter space which closes 
in the thermodynamic limit Q . A scenario of many ACs 
with a broad distribution of widths [134111, as a mani- 
festation of a strong coupling of many energy levels, is 
naturally found in quantum chaotic systems (l3j |. The 
question arises whether the applicability of fidelity can 
be lifted from pure ground-state analysis T6:] to detect 
and characterize ACs in the entire spectrum of a complex 
quantum system. We show in this Rapid Communication 
that the fidelity indeed is a practically useful measure, 
which is local in the spectrum (detecting ACs of a sin- 
gle level) and connects changes in the wave function to 
parametric spectral dynamics. 

Fidelity measure. - Given some parameter depending 
Hamiltonian H(X) = H i + XH 2 , the fidelity @ between 
the n-th eigenstates, denoted by \n), of slightly different 
Hamiltonians H(X) and H(X + SX) is denoted as 

f n (X,SX) = \(n(X)\n(X + SX))\. (1) 

In complex quantum systems with many degrees of free- 



dom, many of the levels of the system are coupled to each 
other leading to ACs in the spectrum of the Hamiltonian 
when the parameter A is changed To simplify the 

discussion, we assume a finite size Hilbert space H, where 
all energy levels are never exactly degenerate. To detect 
and characterize an AC for a given quantum level n we 
study the fidelity change Q 
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which measures the change of the state |n). For 5X <C 1, 
it is independent of SX, i.e. S n (X, SX) « S n (X), and van- 
ishingly small everywhere except in the vicinity of an 
AC. This fidelity measure also has the advantage of be- 
ing applicable locally in the spectrum, where one follows 
a certain state \n(X)) and its neighbors over a range of 
parameter values A to study the ACs they encounter. 
In addition, it is well-suited for numerical computations, 
since A is the only relevant parameter as long as SX is suf- 
ficiently small. The different limit of large SX and hence 
the coupling over a broad energy band was the focus of 
a recent work using another generalized fidelity [l7j . In 
contrast, our interest here is the detection and character- 
ization of ACs as local couplings in energy space. 
Let us first discuss two example systems. An isolated AC 
can locally be described in nearly-degenerate perturba- 
tion theory as an effective two-level system. It is then 
represented by a Hamiltonian H(X) = Xa z + ga x , with 
a real coupling g between the levels [a x and o~ z denote 
Pauli matrices), showing an AC at A = of width c = 2g. 
The normalized eigenstates are easily found fl3l j and the 
fidelity change in the limit <5A 1 can be computed ex- 
actly: 



SiAX) = I 
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(3) 



This is the square of a Lorentzian and differs significantly 
from zero only near the AC at A = 0. This simple ana- 
lytic formula allows us a better understanding of isolated 
ACs, as, for example, the peak width is easily computed 
as cr FWHM = 2g^V2- 1. On the other hand an AC 
can be characterized by the ratio between the local en- 
ergy level curvature and the distance between the two 
repelling energy levels. We call the absolute value of this 
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ratio renormalized curvature C n (A) and find 
1 d 2 E ± {\) 



C±(A) 



A(A) 5A 2 



45±(A) 



(4) 



for the two-level system. For higher-dimensional systems 
we expand the wave function |n(A + 8X)) in second order 
in 5X and find 
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where we reduced the sum near an isolated AC to the 
nearest neighboring level n' . Similarly, one obtains for 
the renormalized curvature 11811 



C„(A) 



2 K"*(A)|ff 2 |n(A))| a 



A(A) 



E n — E v 



J(n'(X)\H 2 \n(X))\ 2 

2 [E n -E n , ? =4 ^ (A) - (5) 



The relation C n « 45„ thus holds as long as the effect of 
other levels can be neglected close to a single AC. 
For very dense spectra and coupling between many levels, 
different ACs may lie close to each other. Three cross- 
ing levels can be generated, e.g., by the following real 
symmetric Hamiltonian 



H(X) 




(6) 



which generalizes the above 2 x 2- model and contains al- 
ready features of more complex systems. Fig. [1] shows 
that the fidelity change, defined in Eq. ([2]), is able to de- 
tect and to distinguish two nearby ACs in this system. 
Furthermore it resembles specific features of an AC in the 
shape of its peak, i.e., depending on the coupling <?, 5n(A) 
shows a narrow peak of height S(X = 0) = l/(8g 2 ). We 
see already in this simple example that the renormalized 
curvature captures the form of the fidelity change 5„(A) 
close to an AC, with deviations arising from the admix- 
ture of a further level, which first and foremost effects 
the local curvature, i.e., the numerator in Eq. (T5|). But 
it also demonstrates that the fidelity change 5(A) itself 
is still effective in detecting and characterizing the ACs. 
Application to complex systems. - A highly dense spec- 
trum with many and possibly overlapping avoided cross- 
ings is encountered in quantum chaotic systems as de- 
scribed by Random Matrix Theory (RMT) [l3[ ■ A prime 
example having such a dense complex spectrum is the 
combination of two random matrices drawn from the 
Gaussian orthogonal ensemble (GOE) [l3| 

H(X) = cos(A)ifi + sin(A)if 2 , H u H 2 e GOE. (7) 

The distribution of minimal distances c at the ACs (nor- 
malized to unit mean) is then given by a Gaussian distri- 
bution P(c) = (2/n) exp [-c 2 /tt] [14]. Using our fidelity 




FIG. 1: (color online) (a) Energy spectrum of Eq. ((6| for 
a — 0, b = 2, c = 3. All levels are coupled and the spec- 
trum shows two close ACs; (b) fidelity change S„(X) and (c) 
renormalized curvature C„(A) for the energy levels of (a), d) 
Energy spectrum for a = 1, b = 2, c = 3. All three levels are 
now directly coupled and the spectrum shows two close ACs. 
(e) S„(X) and (f) C„(X) for the energy levels of (d). 



measure, we can directly detect the avoided crossings in 
this system and estimate also their widths. In the vicin- 
ity of a local maximum, the 5-function has a Lorentzian 
shape as in Eq. ([3]) even in very dense quantum chaotic 
spectra. Under this assumption, we can thus extract the 
width of the AC as c = 2g = \/25 max , c.f. Eq. (0, from 
the local maximum 5 ma x- Averaging over many ACs, 
the fidelity allows the verification of the RMT prediction 
with high accuracy. This is demonstrated in Fig. [2] for 
large random matrices. 

To further exemplify the value of our fidelity measure, 
we apply it to a one-dimensional Bose-Hubbard Hamil- 
tonian with additional Stark force [H Q . This example 
of a many-body Wannier-Stark system can be realized 
with ultracold atoms in optical lattices and the relevant 
parameters may be changed using well-known experimen- 
tal techniques [1| . This model describes N particles on L 
lattice sites, with hopping between adjacent sites and a 
local on-site interaction. As exemplified in [H, H[ , a gauge 
transformation into the force accelerated frame of refer- 
ence turns a constant Stark force into a time-dependent 
phase exp(±iFt) with periodicity Tb = 2ir/F (the Bloch 
period). The corresponding Hamiltonian reads 



l+l 
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u L 
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■1), (8) 



where a j (a/) creates (annihilates) a boson at site I 
and ni — alai is the number of bosons at site I. The 
parameter J is the hopping matrix element, U the 
interaction energy for two atoms occupying the same 
site, and F the Stark force. Periodic boundary con- 
ditions are imposed for H(t), such that the Hamilto- 
nian and the one-period Floquet operator Uf(Tb) = 
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FIG. 2: (color online) Cumulative distribution of ACs deter- 
mined from the fidelity change maxima for the RMT model 
of Eq. © with dimi/i, 2 = 1024 and A G [0,7r[ showing ca. 
30,000 ACs: the numerical distribution (solid line) in excel- 
lent agreement with the RMT prediction CDF(c) = erf(c/-y/7r) 
(dashed line). Inset: Distribution of widths of the ACs P(c) 
(histogram) and the RMT prediction (dashed line). 



Texp ^— i Jq B H(t)dtj (where T denotes time-ordering) 

decompose into a sum of operators for specific quasimo- 
menta k ||. In the following, we use F as a control 
parameter. For J U <C F the quasienergy spectrum 
(eigenphases of Uf(Tb)) is dominated by the force F and 
the system is regular. Decreasing the force to J ps U > F 
the quasienergy spectrum reorders and the coupling be- 
tween the levels becomes more important. For fillings 
of order unity, e.g. N/L ps 1, the system is quantum 
chaotic in this regime and the spectrum obeys Wigner- 
Dyson statistics 0, H| ■ As F is varied one observes an in- 
creasing number of ACs as the spectrum is changing and 
additionally many broad ACs once the quantum chaotic 
region is reached. 

To illustrate the crossover between regions with few 
and many ACs, we study the density of ACs as detected 
by the fidelity change S n , when changing the system pa- 
rameter A. In a histogram, the density Pac(A) is denned 
via pAc(A) • dX = A f Ac(A)/dim'H, comparing the num- 
ber of ACs ./Vac (A) in the interval [A, A + dX] to the total 
number of energy levels dim'H. This is shown in the main 
part of Fig. [3] where we observe no ACs at large F, i.e., 
small values of 1/F, and an increasing number of ACs 
for larger values of 1/F that saturates around 1/F ps 20. 
The mentioned transition between regular and chaotic 
spectral properties for J ps U and F <^ J and approxi- 
mately integer filling in the tilted system can be visual- 
ized by comparing the actual level spacing distribution to 
a Wigner-Dyson distribution using a standard statistical 
X 2 test @. This is displayed in the inset of Fig. [3] along 
with the density of ACs in Fig. [3] The fidelity change 
S(l/F) detects ACs and marks the crossover from regular 
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FIG. 3: (color online) Density of ACs pac in the quasienergy 
spectrum of the Floquet operator of our Bose- Hubbard model 
for varying F and fixed J = 0.038, U = 0.032, N = L = 6. 
The number of ACs as detected by the fidelity change in- 
creases with 1/F and saturates around 1/F w 20 to an aver- 
age value which is shown by the dashed line. Inset: Magnifi- 
cation of the region marked by the box on logarithmic scale 
with a comparison to a v 2 test (with small values for good 
Wigner-Dyson statistics [g|). 



to chaotic dynamics, showing the same qualitative behav- 
ior as the spectral statistics: in regions of good Wigner- 
Dyson statistics we find a high density of ACs compared 
to a smaller number of ACs in the regular regime. The 
crossover beginning for \og(l/F) ps 2, where the density 
of ACs rises above unity, i.e., on average each energy level 
undergoes more than one AC in the unit interval. The 
transition is complete for log(l/F) ss 3 where the \ 2 test 
saturates around a low value. 

A clear advantage of the fidelity change compared to 
spectral statistics is that it can be applied locally in the 
spectrum. This means that, if one is interested only in 
local spectral properties of a system, it is sufficient to 
follow a small number of levels to characterize the sys- 
tem's behavior. On the other hand, it means if we follow 
all energy levels over a certain parameter range, we can 
resolve remarkable details in the full spectrum. This can 
be done easily by using the fidelity change in order to 
detect and characterize ACs. With this method we are, 
e.g., able to detect a small number of regular states [l9| 
traversing the chaotic sea of energy levels in the chaotic 
regime of the tilted Bose-Hubbard model. In this case 
the distribution of widths of ACs is a mixture of regular 
and quantum chaotic distributions: 



P(c) = (1 - 7 )<5(c) + 
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exp 



7 2 c 2 



(9) 



with a chaotic part of weight < 7 < 1 [20j. A finite 
regular component makes itself visible as a strong en- 
hancement of P(c) close to zero, c.f., the inset of Fig. @] 
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FIG. 4: (color online) Cumulative distribution of ACs de- 
termined from the fidelity change maxima for the system of 
Eq. ((8|. Shown are the numerical distribution (solid line), 
the best fit for a mixed RMT spectrum (thick dashed line, 
chaotic part 7 « 0.94), and the RMT prediction for a purely 
chaotic spectrum (thin dashed-dotted line). Parameters: 
N = 6, L = 7, J = 0.038, U = 0.032, F = 1/39 . . . 1/35. Inset: 
Distribution of widths of ACs for the same model (histogram) 
and Eq. ([9]) with 7 w 0.94 (dashed line). The enhancement 
peak close c = arises from the regular "solitonic" states [lj| 
in the spectrum. 

We are also able to estimate the size of this compo- 
nent by analyzing the cumulative distribution function 

CDF(c) = 1 — 7 + 7 erf (-^=^ ■ The result is shown in the 

main part of Fig. 01 where we plot the numerically ob- 



tained distribution and the best x^-fA including a finite 
regular component. We obtain a chaotic part of 7 ~ 0.94, 
corresponding to ca. 6% of regular levels, in good agree- 
ment with counting 7 regular levels out of 132 by direct 
inspection of the spectrum. Except for the identifica- 
tion of single regular levels [l9| , this has so far not been 
detected in the tilted Bose-Hubbard model by other sta- 
tistical measures. The reported results are obtained for 
periodic boundary conditions applied to the Hamiltonian 
of Eq. ([8]), but we found a qualitatively similar picture 
for hard- wall boundary conditions, as used in [19] . Our 
results underline the value of fidelity as a measure for 
detecting ACs with high resolution in complex energy 
spectra. 

Summary. - We showed that quantum fidelity is perfectly 
suited to detect and characterize ACs in the energy spec- 
trum. It therefore connects information about the wave 
function of a system with its spectrum, without direct 
reference to the energy levels by using only the overlap of 
wave functions [2l[ . This has been exemplified for simple 
models and for complex quantum systems showing many 
ACs. The fidelity of Eq. ([2]) can be applied to charac- 
terize the system's behavior locally in energy space, i.e., 
parts of spectrum only, as is the typical situation in ex- 
periments where the complete spectrum is not accessible 
in most cases. The fidelity, therefore, proves very useful 
to study many-body systems, also beyond their ground- 
state properties Q. 
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